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ABSTRACT 

We present an optical analysis of a sample of 11 clusters built from the EXCPRES sample of X-ray selected clusters at intermediate 
redshift (z ~ 0.5). With a careful selection of the background galaxies we provide the mass maps reconstructed from the weak lensing 
by the clusters. We compare them with the light distribution traced by the early-type galaxies selected along the red sequence for 
each cluster. The strong correlations between dark matter and galaxy distributions are confirmed, although some discrepancies arise, 
mostly for merging or perturbed clusters. The average M/L ratio of the clusters is found to be: MIL,- = 160 ± 60 in solar units (with 
no evolutionary correction), in excellent agreement with similar previous studies. No strong evolutionary effects are identified even if 
the small sample size reduces the significance of the result. We also provide a individual analysis of each cluster in the sample with 
a comparison between the dark matter, the galaxies and the gas distributions. Some of the clusters are studied for the first time in the 
optical. 

Key words. Gravitational lensing: weak - X-rays: galaxies: clusters - Cosmology: observations - Cosmology: dark matter - Galaxies 
: clusters : general - Galaxies : clusters : individual 


1. Introduction 


Although the origin and evolution of linear-scale clustering is 
well described by the concordance model ( |Spergel et al.|2007| , 
gravitational clustering of matter on smaller scales (galaxy clus¬ 
ters and groups) belongs to a non-linear regime of structure for¬ 
mation. This regime is more difficult to understand and to sim¬ 
ulate because its evolution must include the role of baryons, 
driven by complex physics. Clusters of galaxies which are 
the most massive gravitationally bounded structures have been 
widely used over the past years to probe the cosmic evolution of 
the large scale structures in the Universe ( |Voit|2005 Allen et al. 
|201 l| l. In the standard model of structure formation driven by 
gravitation only, clusters form a self-similar population charac¬ 
terized only by their mass and redshift. Including baryon physics 
introduces some distortions in the scaling relations between the 


* Based on observations obtained with MegaPrime/MegaCam, a joint 
project of CFHT and CEA/DAPNIA, at the Canada-Erance-Hawaii 
Telescope (CFHT) which is operated by the National Research Council 
(NRC) of Canada, the Institut National des Science de I’Univers of the 
Centre National de la Recherche Scientifique (CNRS) of France, and 
the University of Hawaii. This research also used the facilities of the 
Canadian Astronomy Data Centre operated by the National Research 
Council of Canada with the support of the Canadian Space Agency. 
Also based on observations obtained with XMM-Newton, an ESA sci¬ 
ence mission with instruments and contributions directly funded by 
ESA Member States and NASA. 


mass and other physical quantities like temperature, X-ray or 
optical luminosity, etc ... ( Kaiser| 1986[ Giodini et al.| 20T3] l. 
Most research works have recently focused on the relationship 
between the dominant dark matter and the baryonic matter that 
forms gas and stars ( Lin et al.|2003[ Giodini et al.|200^ . Both 
the mass-to-light (M/L) ratio of structures and the halo occu¬ 
pation number (HON, or the number of satellite galaxies per 
halo) correspond to observables easy to compare to predictions 
from numerical simulations (Cooray & Sheth|2002t Tinker et al. 


20051. They are both representative of the way stellar forma¬ 


tion occurred in the early stages of halo formation (Marinoni & 
Hudson 2002[ Borgani & Kravtsov 201 l|l. Recent progress on 


numerical simulations ( Murante et al.|2007] [Conroy et al.|2007 
Aghanim et al. 2009|l have also stressed the role of hierarchi 


cal building of structures in emdching the intra-cluster medium 
(ICM) with stars in a consistent way with the observed amount 
of ICM globular clusters and ICM light. This ICM light, al¬ 
though hardly detectable, can be considered as the extension of 
the diffuse envelope often seen in the central galaxy in rich clus¬ 
ters of galaxies. It is an important component, although not the 
only one, which explains the formation of the Brightest Cluster 
Galaxies (hereafter BCG) in the centre of clusters of galaxies 
( |Dubinski|1998[[Presotto et al. 12014) 1. 


To quantify these processes and to compare them with those 
included in numerical simulations, it is of prime importance to 
get reliable masses and mass distribution in clusters of galax- 
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ies. But getting accurate mass estimates is a difficult task and 
large uncertainties reduce the validity of the relation between the 
light and total mass distribution (Vale & Ostriker 2004[ Rozo 
|et al.|[2OT4) l. The determination of the mass distribution using 
the weak gravitational lensing of the background galaxies by 
clusters of galaxies is a powerful approach to address this ques¬ 
tion ( [Schneider et al.|[2006[ [Hoekstra et al.||T013j l. Lensing is 
able to trace directly the dark matter component in rich clus¬ 
ters of galaxies down to massive groups: see a few examples in 
[Broadhurst et al. ( 2005| l ;|Gavazzi & Soucail ( |2007]l ; |Bradac et al. 


( 2008| l; Umetsu et al.| ( |2012| l; Zitrin et al. ( |2013| l; Gastaldello 
et al. (2014|l. Recent analysis of cluster samples have shown 


strong improvements in the accuracy of the mass measurements 


(jvon der Linden et al.|2014[ Applegate et al.|2014| Kettula et al. 


2014| l, thanks to the powerful capacities of ground based wide 
field imaging. Finally, spectacular results were obtained in the 
case of merging clusters like the “Bullet” cluster for which the 
dark matter distribution traces closely the galaxy distribution 
while the intra-cluster gas traces by X-ray emission is rather 
uncorrelated from the non-collisional components (jClowe et al. 
|2006l l. 


Another way to characterize the relation between the mass 
distribution and the stellar light is to compare the M/L ra¬ 
tio with cluster properties. In particular, recent results obtained 
with weak lensing masses seem to show a slight scaling depen¬ 
dence of the M/L ratio with mass: this is demonstrated from the 
MaxBCG sample built from the SDSS, with structures ranging 
from small groups to massive clusters (Sheldon et al.|20()9ll and 


also from samples of clusters of galaxies ( jMuzzin et al.pOOT 
Popesso et al.|2007 Bardeau et al.|2007 1. All these clusters are 


at low redshift (< 0.2 typically) because there are severe obser¬ 
vational limits at larger redshifts. 


The purpose of this work is to present a detailed view of the 
relations between dark matter and stellar light from a sample of 
11 clusters at intermediate redshift (~ 0.5). This sample is part 
of the EXCPRES sample of cluster, built as an un-biased sam¬ 
ple of clusters of luminous X-ray clusters at redshift around 0.5, 
covering a wide range of dynamical mass and X-ray tempera¬ 
ture. The whole sample of 29 clusters was observed in X-ray 
with XMM- Newton (Arnaud et al. in preparation) to test the 
evolution of clusters properties with redshift. The 11 brightest 
ones were observed at CFHT for optical follow-up. A weak lens¬ 
ing analysis was proposed to provide a mass estimate for these 
clusters. Its practical implementation, as well as the global mass 
analysis were presented in Foex et al. (2012 hereafter Paper 1). 
In the present paper, we focus on the comparison between the 
optical properties of the clusters and their mass distribution and 
we present the characteristics of each individual cluster. The pa¬ 
per is organized as follow : Section 2 presents the data used in 
the analysis and the selections of the different catalogs. Section 
3 presents the global optical properties of the clusters while 
Section 4 is dedicated to the dark matter bi-dimensional distri¬ 
bution from the weak lensing maps. In Section 5 we discuss the 
properties of the sample and the links between the stellar light 
distribution and the total mass. Conclusions are given in Section 
6. The individual properties of the 11 clusters of the sample are 
detailed in the Appendix. 


Throughout this paper, we use a standard A-CDM cosmol¬ 
ogy with Qm = 0.3, Ga = 0.7 and a Hubble constant Hq - 70 
km s“' Mpc^' or h = //q/IOO = 0.7. 


2. Observational data 

2.1. Observations and data reduction 


Data were obtained for the whole cluster set with the MegaPrime 
instrument at the Canada-France-Hawaii Telescope, during the 3 
observing periods in 2006 and 2007 (RunIDs: 06AF26, 06BF26 
and 07AF8; PI: G. Soucail). The camera MegaCam is a wide 
field CCD mosaic covering 1 square degree, with a pixel size 
of 0.186". Multi-color imaging was obtained with the 4 photo¬ 
metric filters g', r', i' and z' with integration times of 1600 sec., 
7200 sec., 1200 sec. and 1800 sec. respectively. In practice, some 
clusters were observed with slightly longer integration times due 
to the re-observation of some images in poorer observing condi¬ 
tions. All data obtained in r' were done in good seeing conditions 
(IQ less than 0.8") and during photometric periods. The integra¬ 
tion time in r' was defined to obtain a limiting magnitude for 
weak lensing studies / ^ 26: we consider that at this magnitude 
limit, the bulk of the background sources is at redshift higher 
than 1 and that the lensing strength of the clusters is maximum. 
For the 3 other colors, the strategy was dehned to detect cluster 
galaxies up to m* + 4 (with m* ^ 20 at z = 0.5) in reasonable 
seeing conditions (IQ < 1"). The summary of the observations 
is presented in Table 

Data reduction was done in a standard way for large CCD 
mosaics. After the on-line preprocessing done at CFHT (cor¬ 
rection of instrumental pixel-to-pixel effects) with the Elixir 
pipelin^ the global processing was done either by the Terapix 
team or by the authors, using the Terapix tool^locally. A global 
astrometric solution was found with Scamp using the star refer¬ 
ences of the USNO-B1 catalog ( [Monet et al.|2003| l, as well as the 
photometric alignment of the different images. The final stack¬ 
ing of each image set was done with Swarp, producing a single 
wide field image and its associated weight map. A;(f^-image was 
built from g',r',i' images and was used for objects detection. 
Note that efficient flat-fielding in the far-red was difficult so the 
z'-images were not included in the;\f^-images. 

The cluster RXJ1347.5-1145 was retrieved directly from the 
CFHT-CADC archive^ it was observed in g' and r' (PI: H. 
Hoekstra, runID: 05AC10) and extensively analyzed previously 
in weak and strong lensing ([Bradac et al.[2008| Halkola et al. 


[2008[ ). Note also that the field of view of ELXJ2228.5-1-2036, be¬ 
ing at low galactic latitude, is highly contaminated by bright 
stars. It is necessary to mask large areas around all the stars and 
this prevents an efficient weak lensing analysis. Results obtained 
with this cluster will have to be taken with caution. 


2.2. Muiti-coior photometry 

Photometric catalogs were built with Sextractor (Bertin & 
[Arnouts[[1996[ l in dual mode, the detection of the objects be¬ 
ing made on the -image for each cluster. For the magnitudes 
of the objects, we used the MAG^UTO parameters while for 
color indices, we used the MAG_APER magnitudes measured 
in an constant circular aperture of 3" in diameter. We also cor¬ 
rected uncertainties in the zero point calibration, a critical step 
for further estimates of photometric redshifts. We used the color- 
color distributions of the stars detected in each field and we 
compared these distributions with the expected ones computed 

’ http://WWW.cfht.hawaii.edu/Instruments/Elixir/ 

^ http://www.terapix.fr/ 

^ The Canadian Astronomy Data Centre is operated hy the National 
Research Council of Canada with the support of the Canadian Space 
Agency. 
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Table 1. General properties of the clusters and summary of the observations done in the r' band, i.e. the data used for the weak 
lensing analysis. The last columns give the mean galaxy number densities before and after removal of foreground and cluster 
contamination (in units of arcmin^^). 


Cluster 

RA 

(J2000) 

Dec 

(J2000) 

Redshift 

Exp. time 
(sec.) 

Seeing 

(") 

Completeness 
in r' (50%) 

Galaxy 

density 

Background 
galaxy density 

MS 0015.9±1609 

00'T8'"33.26' 

±16°26'12.9" 

0.541 

5600 

0.82 

24.50 

26.0 

18.0 

MS 0451.6-0305 

04''54'’T0.85' 

-03°00'57.0" 

0.537 

7200 

0.77 

24.75 

30.6 

23.5 

RXC J0856.1±3756 

08''56'’T2.69' 

±37°56'15.0" 

0.411 

7200 

0.66 

24.90 

32.7 

23.3 

RX J0943.0±4659 

09'’42'”56.60' 

±46°59'22.0" 

0.407 

7200 

0.87 

24.60 

23.7 

16.5 

RXC J1003.0±3254 

10''03'”04.62' 

±32°53'40.6" 

0.416 

7200 

0.79 

24.55 

26.2 

19.2 

RX J1120.1±4318 

ll''20'”07.47' 

±43°18'06.0" 

0.612 

7200 

0.60 

24.85 

29.6 

23.5 

RXC J 1206.2-0848 

12''06"T2.13' 

-08°48'03.6" 

0.441 

7200 

0.85 

24.90 

32.0 

26.5 

MS 1241.5±1710 

12''44'”01.46' 

±16°53'43.9" 

0.549 

7200 

0.72 

24.85 

33.0 

16.3 

RX J1347.5-1145 

13''47'"32.00' 

-ir45'42.0" 

0.451 

7200 

0.77 

24.95 

29.6 

25.8 

MS 1621.5±2640 

16'’23'"35.16' 

±26°34'28.2" 

0.426 

7200 

0.60 

25.05 

37.1 

28.9 

RX J2228.5±2036 

22''28'"33.73' 

±20°37'15.9" 

0.412 

7200 

0.69 

24.85 

33.5 

25.6 


by convolving a well calibrated spectral stellar library (Pickles 
1998|l with the biter and instrumental transmissions. The posi 


tion of the “knee” seen in the stellar color-color diagrams was 
also matched to the observed colors. We bnally considered the 
r'-band photometry as a reference and computed the best cor¬ 
rections to apply for the other biters, using a minimization 
between both distributions. 

Separation between stars and galaxies were obtained by fol¬ 
lowing the methodology developed by Bardeau et al. ( |20051 l and 
all stellar-like objects were removed from the catalogs. The 50% 
completeness limit of the galaxy catalogs in the r'-band are given 
in Table [T] 


2.3. Photometric redshifts and cluster member selection 

To identify cluster members in the photometric catalogs, we used 
an updated version of the public code HyperZ|^( Version 11, June 
2009) and computed photometric redshifts. HyperZ is based on 
a template btting method ( Bolzonella et al.||2000[ |Pell6 et ah 
|2009| l: the measured specttal energy distribution (SED) is btted 
with a library of templates built from different spectral types, star 
formation histories and redshifts. In the present case, photomet¬ 
ric redshifts were computed in the redshift range [0,4]. We did 
not try to compute Bayesian redshifts with a prior on luminos¬ 
ity but we used a simple cut in the permitted range of absolute 
magnitudes; -25 < M < -14. 

HyperZ can be very efficient with a sufficient number of pho¬ 
tometric bands but in this work we only have magnitudes in 4 
biters, with different limiting magnitudes so it is quite challeng¬ 
ing to perform good photometric redshifts in all redshift ranges. 
So we did not try to assign a photometric redshift to each indi¬ 
vidual galaxy and we focused mainly on the selection of cluster 
members and the detection of the cluster over-density. To have 
a more quantitative estimate of the reliability of these photo- 
mettic redshifts, we simulated a photometric catalog with sim¬ 
ilar properties as our present observations. For simplicity, we 
used a bat redshift distribution, which is sufficient to test our 
redshift ranges of interest. We generated simulated magnitudes 
for 100,000 galaxies in the 4 MegaCam biters adding noise ac¬ 
cording to the average signal-to-noise ratio observed in our data. 
Then we ran HyperZ on this simulated catalog and compared the 
photometric redshifts to the expected values (Fig.[T]). Because of 
several color-color degeneracies, the reliability of HyperZ is not 


http://webast.ast.obs-mip.fr/hyperz/ 



Fig. 1. Photometric redshifts estimated by HyperZ for a simu¬ 
lated bat distribution of galaxies. The color scale shows the den¬ 
sity of points, from black to yellow. 


constant across the whole redshift range and many high redshifts 
galaxies with Ztme >1.5 have an under-estimated photometric 
redshift Zphot ~ 0.5. But for galaxies with 0.4 < Ztme < 0.6, the 
results are quite satisfactory; most of the galaxies have a “cor¬ 
rect” photometric redshift (±0.1) and only a small fraction of 
them have “catastrophic” values. Therefore, we consider that the 
photometric redshifts estimated by HyperZ can be used safely to 
pre-select cluster members. A consistency check was done by 
looking at the clusters over-density of galaxies; we determined 
the redshift distribution far from the cluster centre and subtracted 
it from the one obtained in a central region covering the same 
area. The resulting redshift distribution shows clearly a peak lo¬ 
cated close to the cluster spectroscopic redshift (Figure]^. 
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Fig. 2. Over-densities in the photometric redshifts distribution for each cluster. In each panel the vertical line shows the spectroscopic 
redshift of the cluster. The plot represents the redshift distribution of the central area defined as R < 5' minus the distribution of an 
annulus of same area starting at /? = 10'. 


2.4. The cluster color-magnitude diagram 


Early-type galaxies form a homogeneous population whose 
spectral energy distribution is dominated by red and old stars. 
For a given redshift these galaxies are distributed along a well 
defined “red sequence” in a color-magnitude diagram and this 
sequence extends over several magnitudes with a small scatter 
(smaller than 0.1 mag typically). This characteristics has long 
been used as a powerful tool to detect clusters of galaxies in large 
photometric surveys ( Gladders & Yee|2000 1 . In the present study 
we used the red sequence of the galaxies to clean the lensing 
catalogs from cluster members. In order to identify this red se¬ 
quence we followed the method described by Stott et al. ( 2009| l: 
we first selected a sub-sample of objects located in the central 
area of the images (R < 200" from the cluster centre) and we 
only kept galaxies with a photometric redshift compatible with 
the cluster redshift. We then performed a linear fit of the red 
sequence with a 3cr clipping to determine analytically the color- 
magnitude relation in the (r',r' - i') diagram (Fig. [^. On av¬ 
erage, the dispersion cr around the red sequence is 0.07 + 0.02. 
Finally we excluded all the galaxies in the whole catalog along 
this relation and within +3cr of the Gaussian fit. We also applied 
a magnitude cut 18 < < 23 because fainter galaxies are no 

longer dominated by cluster members and it is of prime impor¬ 
tance to keep a background density as large as possible for the 
weak lensing analysis. 

In summary, the background galaxies catalogs were built 
for the present work with the following rules: a magnitude cut 


22 < r' < 26 and a color cut outside the red sequence +3(T, up 
to r' = 23. Thus, most of the cluster members were removed. 
We checked that the galaxies density profile of the remaining 
galaxies is flat, except very near the cluster centre where some 
residual contamination remains (see Paper 1). However, this has 
no strong impact of the global morphology of the clusters that 
is described in this paper. The values of the average background 
galaxies density are given in Table before and after cleaning 
the catalogs. 


3. Stellar light distribution 

3.1. Selection of cluster galaxies and global cluster 
properties 


As stated in Section 2, we specifically built galaxy cluster cata¬ 
logs by selecting galaxies for which the color r' - i' falls within 
+3cr of the cluster red sequence. We also limited the sample 
to galaxies brighter than 0.4L*, i.e. m* H- 1 to avoid too much 
contamination in the faint magnitude bins. The m* magnitude 
was computed for each cl uster, assuming an abs olute magnitude 
M* - 5 log h - -20.44 (Blanton et al. 2003 i and adding the 
appropriate distance modulus, the k-correction for an early-type 
galaxy for each cluster and the galactic reddening correction. 
Adding an evolution correction in the luminosity of the ellipti¬ 
cal galaxies would amount to ~ 0.65 mag. This would decrease 
the intrinsic luminosity at a given apparent magnitude by a fac¬ 
tor 1.8 and it would also change the magnitude cut in the galaxy 
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Table 2. Global properties of the cluster sample: the radius R 200 and the 2D projected mass are derived from the weak lensing 
analysis. The optical luminosity L 200 and the galaxy number N 200 are measured inside the radius R 200 and are corrected from the 
background contamination. The total luminosity has been corrected from the incompleteness at faint magnitudes. For comparison 
we also give the number of cluster members inside a physical radius of 1 Mpc, N (IMpc). The global M/L ratio is given in solar 
units in the r'-band. 


Cluster 

R 200 

(h-^ Mpc) 

^^^200 

T 200 

(lO'^h-^LQ) 

N 200 

N(lMpc) 

M/L 

(/i7oMo/LO) 

MS 0015.9+1609 

2.33 

±0.13 

3.27 ± 0.58 

17.2 ± 1.6 

135 ± 12 

69 ±8 

190 ±51 

MS 0451.6-0305 

1.92 

±0.12 

1.84 ±0.35 

10.8 ± 1.2 

81 ±9 

48 ±7 

170 ±52 

RXCJ0856.1+3756 

1.65 

±0.10 

1.0 ±0.175 

5.1 ±0.7 

48 ±7 

28 ±5 

197 ±63 

RX J0943.0+4659 

1.77 

±0.80 

1.18 ±0.25 

13.6 ± 1.4 

100 ± 11 

64 ±8 

87 ±28 

RXC J1003.0+3254 

1.69 

±0.11 

0.94 ± 0.20 

4.8 ± 0.9 

28 ±5 

29 ±5 

198 ±81 

RXJl 120.1+4318 

1.46 

±0.14 

0.75 ± 0.26 

8.2 ± 1.2 

52 ±7 

46 ±7 

91 ±45 

RXC J1206.2-0848 

2.03 

±0.10 

1.93 ±0.29 

17.4 ± 1.7 

111 ± 11 

54 ±7 

111 ± 27 

MS 1241.5+1710 

1.88 

±0.13 

1.78 ±0.38 

8.0± 1.1 

55 ±8 

30 ±6 

224 ± 79 

RX J1347.5-1145 

2.40 

±0.10 

3.27 ± 0.40 

13.9 ± 1.4 

99 ± 10 

42 ±7 

235 ± 53 

MS 1621.5+2640 

1.90 

±0.11 

1.57 ±0.26 

7.9 ± 1.0 

63 ±8 

37 ±6 

199 ±59 

RX J2228.5+2036 

1.68 

±0.12 

1.09 ±0.24 

16.0 ± 1.9 

117± 13 

58 ±8 

68 ±23 


m 



Fig. 3. Color-magnitude diagram in the held of 
MS 1621.5+2640. The black points are the galaxies located 
at r < 500" from the cluster centre. The red dots are those with 
r < 200" and a Zphoi compatible with the cluster redshift. The 
straight line is the best ht to the red sequence and the dashed 
lines are the 3cr limits of the red sequence. At the cluster redshift 
{z = 0.43), the expected color, computed with the synthetic 
evolutionary code of Bruzual & Chariot ( 2003|l is r' - i' - 0.77. 


catalogues and reduce the number of cluster galaxies. All in all 
the expected change in the total luminosity of the clusters is a 
factor of 2 to 2.2. But it is highly uncertain because it strongly 
depends on the galaxy evolutionary scheme used in the mod¬ 
elling of the evolution correction and it is generally not included 
in the global studies of clusters. To remain consistent with pre¬ 
vious works (see |Popesso et al. (20041; Bardeau et al. (20071 for 
example), we chose not to include it in the present work. 

The total luminosity of the clusters was measured by sum¬ 
ming the luminosity of all galaxies in the red sequence interval 
and located inside the radius /?200 estimated with the weak lens¬ 
ing analysis (Paper 1). A background correction to the total lu¬ 
minosity was included by removing an average luminosity mea¬ 
sured in an annulus dehned by 2 /?200 < r < 3 /?200 for each cluster 


and scaled with the adequate area. This correction is rather small 
because the total cluster luminosity is dominated by bright early- 
type galaxies. Finally a correction for the magnitude cut of the 
galaxy catalogs is added (Popesso et al. 2004| l. It is calculated as 
the integral of a Schechter function up to 0.4L* and with a slope 
a - 1.25 ( [Blanton et al.|20()^ . A factor 1.6 is then included to 
get the total luminosity of the clusters (Table |^. We also com¬ 
puted the cluster optical richness A 200 which we dehne as the 
number of cluster galaxies within f ?200 after correction from the 
background ([Hansen et al.|2005|l. All values are given in Table 

m -— 


3.2. Morphology of the light distribution 

There are many possibilities to map the light distribution in a 
galaxy cluster. They depend on the choice of the input catalogs 


and the method used to derive a correct mapping (see Okabe 
[et al.[ ( [2010[ l for example). In the present case we generated the 
light map trying to take into account more accurately the elliptic- 
ity and the orientation of each individual galaxy for the building 
of the cluster light distribution. In practice we selected cluster 
galaxies with magnitudes 18 < < 23 and colors along the red 

sequence (±3cr). An artihcial image was created for each with 
the package Artdata in IRAF|^ with parameters extracted from 
the photometric catalogs of the cluster galaxies (position, mag¬ 
nitude, ellipticity e and position angle PA). The resulting image 
was then smoothed with a Gaussian kernel to generate the cluster 
light map. We fixed the FWHM of the kernel to 80", a size two 
times smaller than the smoothing scale of the dark matter mass 
map. The maps are displayed for each cluster together with the 
mass distributions (Figs. A.l to A.l 1[). The ellipticity e - l-bja 


and PA of the cluster light were measured on a 2D fit of the 
isophotes with Ellipse. For each cluster we favored the large 
scale morphology and we htted isocontours which in some cases 
encompass several clumps, especially in clusters with complex 
structure. In most cases this corresponds to a radius 100 to 150", 
i.e. up to 1 Mpc at the cluster redshift (Table [^. Error bars on 

^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 


5 


















































G. Soucail et al.: The matter distribution in z ~ 0.5 redshift clusters of galaxies 


Table 3. Morphological properties of the luminous component of the clusters. The ellipticity e — 1 -bja and the position angle PA 
are given for the clusters luminosity density map as well as the semi-major axis of the measured ellipse. For the Brightest Cluster 
Galaxy (BCG), only e and PA are given. PA are in degrees, counted counter-clockwise with respect to the North-South axis. 


Cluster 

6 

Cluster light 

PA (deg.) radius (kpc and ") 

BCG 

6 

PA (deg.) 

MS 0015.9-P1609 

0.20 ±0.10 

36 ±9 

520 (80") 

0.26 ± 0.03 

75 ±8 

MS 0451.6-0305 

0.24 ± 0.02 

150 ±3 

510(89") 

0.28 ± 0.02 

101 ±8 

RXC J0856.1-1-3756 

0.23 ± 0.04 

2±5 

660 (132") 

0.28 ± 0.01 

135 ± 10 

RX J0943.0-P4659 

0.27 ± 0.05 

48 ±2 

470 (95") 

0.28 ±0.01 

11 ± 10 

RXC J1003.0-P3254 

0.53 ± 0.03 

22 ±3 

450 (89") 

0.28 ± 0.02 

-2± 10 

RX J1120.1-P4318 

0.38 ± 0.02 

95 ±4 

720(119") 

0.40 ± 0.04 

110± 19 

RXC J 1206.2-0848 

0.36 ±0.01 

86 ±3 

680(130") 

0.49 ± 0.02 

111 ± 24 

MS 1241.5-P1710 

0.10 ±0.03 

79 ±4 

260 (45") 

0.16 ±0.01 

5± 10 

RX J1347.5-1145 

0.36 ± 0.08 

43 ±3 

430 (80") 

0.31 ±0.03 

0± 10 

MS 1621.5-P2640 

0.27 ± 0.04 

22 ±8 

380 (74") 

0.02 ± 0.02 

93± 14 

RX J2228.5-P2036 

0.30 ±0.10 

69 ±5 

280 (56") 

0.30 ± 0.03 

93± 14 


0.8 


O 0-6 
o 

I 

p. 
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Ellipticity (Optical light) 

Fig. 4. Ellipticity of the global light distribution of the clusters 
versus the ellipticity of the brightest central galaxy. 
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the elliptical parameters are estimated from the change of the 
parameters when the radius varies from 90 to 150 " typically. 

3.3. The brightest cluster galaxy (BCG) 

BCGs are usually located at the very centre of clusters of galax¬ 
ies. Numerous studies emphasized their specific properties com¬ 
pared to lower luminosity cluster me mbers (|Lin & Mohr|20(j^ 
[Smith et al.|2010HHaarsnia et al.|2010[|Ascaso et al |201 1| |: lumi- 
nosity, size and effective radius, star formation history and stel¬ 
lar populations ... Discussions to discriminate between the role 
of internal feedback processes and of environment and merging 
of satellite galaxies in the formation of the BCG are still open. 
They also depend on the scenario of galaxy formation and the 
ACDM paradigm seems to favor the importance of galaxy merg¬ 
ers at the centre of the main halo ( jPe Lucia et al.|2007| l. On the 
contrary recent observations confirm the importance of baryonic 
feedback in the size evolution of the BCGs ( jAscaso et al.|201 f) . 
It is out of the scope of this paper to produce a detailed analy¬ 
sis of the structural parameters of the BCGs in our sample, and 
we simply compared their ellipticity and orientation (measured 
with Sextractor) with the large scale light and dark matter dis¬ 
tributions (Table [^. The link with the dark matter ellipticity is 


not obvious but the correlation between the light distribution and 
the central BCG is quite convincing (Fig|^, except for two out¬ 
liers which have complex sub-structures. Quantitatively we find 
a weighted mean < sugkt - sbcc >- 0.001 + 0.12 and the mean 
orthogonal deviation from the 1;1 line [sugi,, = sbcc) is 0.06. 
This value decreases to 0.03 if we remove the two major outliers, 
well below the uncertainties in the ellipticity measurements, of 
the order of 0.05. This result is not surprising as the alignment of 
the BCG with the distribution of galaxies at large scale was first 
observed by Lambas et al.jl 1988|l and confirmed since then b y 
many studies ( Panko et al.|2009[|Niederste-Ostholt et al.|20T0 l. 

In a second step we separated the clusters in two classes: 
those dominated by a single giant elliptical galaxy, most often 
embedded in an extended envelope (a cD-type galaxy) and those 
for which more than one bright galaxy forms the cluster centre, 
or the brightest cluster member does not outshine other galax¬ 
ies. 5 out of 11 clusters are dominated by a cD galaxy, namely 
RXC J0856.1H-3756, RX J1120.1 -h 4318, RXC J1206.2-0848, 
MS 1241.5-H1710 and RX J2228.5-H2036. Surprisingly they are 
not necessarily the brightest clusters in terms of total stellar lu¬ 
minosity nor the most massive ones, suggesting that the forma¬ 
tion of a giant cluster galaxy is not only related to the initial halo 
conditions but also to the evolution processes in the clusters and 
the merging history of the structures (|Dubinski|1998|). 


4. Bi-dimensional weak lensing analysis and dark 
matter distribution 

4.1. Mass reconstruction 


We refer to Paper 1 for the details of the weak lensing implemen¬ 
tation. In summary, the galaxies shapes are measured with the 
Im2shape software ( Bridle et al.||2002] l. For each object a para¬ 
metric shape model is set up with an ellipse. Im2shape convolves 
this model with the local PSF and subtracts it to the sub-image 
centered on the galaxy. A MCMC minimizer applied on the im¬ 
age residuals provides the intrinsic shape parameters and error 
estimates. The PSF is measured directly on the images by aver¬ 
aging the shapes of the 5 closest stars to each galaxy. All mea¬ 
sures are done on the r' images which were obtained with the 
highest image quality. Only galaxies for which the measured el¬ 
lipticity error is smaller than 0.25 are kept in the working cat¬ 
alogs. The main results regarding the mass measurements and 
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the study of the global properties of the clusters are presented 
in Paper 1, as well as a careful analysis of the sources of error 
in the mass determination. In the present work, we focus on the 
spatial distribution of the dark matter traced by the weak lensing 
map reconstruction. 

We used the software LensEnt2 kindly provided by P. 
Marshall ( [Marshall et al.|[2002] l to build the weak lensing mass 
maps. The method is based on an entropy-regularized maximum- 
likelihood technique. It uses the shape of each background 
galaxy as an individual estimator of the local reduced shear. The 
pixel size on the mass grid is chosen to have more or less 1 
galaxy per pixel leading to a comparable number of data and free 
parameters. To consider the fact that clusters have an extended 
and smooth mass distribution, the code includes a smoothing via 
the size of the Intrinsic Correlation Function (ICF): the physical 
mass map is expressed as a convolution of the “hidden” distri¬ 
bution (the pixels grid) with a broad kernel dehned by the ICF. 
The shape and size of this ICF are the main control parame¬ 
ters of LensFnt2 and reflect the spatial resolution of the recon¬ 
structed mass map. For simplicity, we used only a Gaussian ICF 
with a width of 150" for all clusters. This represents a good 
compromise between smoothness and details in the mass map. 
LensFnt2 also computes error maps which give locally the width 
of the probability function in the mass reconstruction. The aver¬ 
age value of the error map within the central area of the CCD 
image (15' x 15') is a good estimate of the level of uncertainty in 
the mass reconstruction. This is the value which determines the 
signal-to-noise of the detected peaks (Table |^. 

LensFnt2 provides output mass maps in physical units of 
surface mass density (Mq pc^^). This is valid provided that the 
redshift distribution of the background sources is well known. 
In the present case, we worked with source catalogs which were 
cleaned for galaxy cluster members, but still contaminated by 
foreground galaxies. A rough estimate of such contamination 
comes from the redshift distribution of galaxies in the magnitude 
range selected for our catalog: if we apply the same selection cri¬ 
teria on the deep photometric catalog with photometric redshifts 
built from the CFHTLS-Deep survey ( Coupon et al.||20()9) l, we 
hnd that about 25% of the galaxies are at redshift smaller than 
0.5, i.e. foreground galaxies. This is coherent with the number 
found in deep spectroscopic surveys although at slightly brighter 
magnitudes ( |Le Fevre et al.|2005] l. In our case, the effect of such 
uniform contamination is mostly a dilution of the weak lens¬ 
ing signal, so the output mass densities of LensFnt2 are not 
reliable in their absolute values. Moreover, there remains some 
additional contamination from cluster galaxies in the very cen¬ 
tre of the clusters (see Figure 3 in Paper 1). The main effect is 
to attenuate the peak intensity and to decrease the S/N ratio of 
the cluster component. But we do not expect any signihcant in¬ 
fluence on the shape of the mass reconstruction, provided the 
contaminating cluster members are randomly oriented within the 
cluster. This assumption is valid in our case because the galaxy 
catalogs include cuts in color and magnitude which eliminate 
all bright cluster members. The remaining galaxies are mostly 
blue and/or faint so they are less sensitive to intrinsic alignment 
effects ( [Mandelbaum et aT]|2011[ ). Therefore in the rest of the 
paper we concentrate our work on the 2D mass distribution. The 
mass map reconstructions are displayed for each cluster in the 
Appendix. 

4.2. Ellipticity of the mass distribution 

The ellipticity of the mass distribution traced by the weak lens¬ 
ing mass reconstruction has been the focus of several studies. It 


is a direct evidence of the triaxiality of the cluster halos and is 
expected to be a non-negligible factor in the growth of massive 
halos in the ACDM paradigm ( [Limousin et al. 20131. Oguri et al. 
( [2010| l studied a sample of 25 massive X-ray clusters (mostly in 
the LOCUSS sample) at redshift ~ 0.2. They found an average 
ellipticity (e) - 1 - bja - 0.46 + 0.04. More recently Oguri et al 


( |2012| l conhrmed the trend with another independent sample of 
clusters built from the Sloan Giant Arcs Surveys as part of the 
SDSS ( Hennawi et al.|2008 i. These values of the mean elliptic¬ 
ity are in good agreement with the theoretical predictions based 
on numerical simulations of cluster dark matter halos (Jing & 
Suto 2002|l. They correspond to what is expected for massive 


clusters, contrary to low mass clusters which are expected to be 
more circular. 

We tried to explore this question with the present sample. 
But because our sample is at higher redshift than LOCUSS, sev¬ 
eral difficulties limit the outcomes of the approach. In order to 
get a estimate of the uncertainties of the elliptical parameters 
for each mass reconstruction, we used a “jackknife” resampling 
method to remove 10% of the galaxies in the source catalog, 
and we repeated the process 10 times. The removed galaxies all 
differ from one attempt to the others. 10 new mass maps were 
computed for each cluster with these sub-catalogs, as well as the 
corresponding error maps. The 10 error maps were averaged and 
the average level of this frame gives the Icr level of the mass 
reconstructions. We then fitted each of the 10 mass maps with 
elliptical contours and selected the elliptical parameters of the 
3cr isocontour. An average of the 10 hts gives the hnal values 
for the elliptical parameters (s and PA), as well as their stan¬ 
dard deviation (Tablej^. This process was acceptable except for 
the clusters RX J1120.1+4318 and MS 1241.5+1710 with the 
lowest S/N maps. We restricted their fits to the 2cr and 2.5cr iso¬ 
contours respectively. 

In all cases this corresponds to an isocontour of 100 to 150" 
in radius (or 600 to 900 kpc at redshift 0.5). But because 
of the limited resolution of the mass reconstruction, the effect 
of the central smoothing by the ICF is quite signihcant and in¬ 
duces an attenuation in the measure of the mass ellipticity. In 
practice, an ICF of 150" corresponds to a Gaussian smoothing 
with cr ~ 60". As a test case we simulated a set of mass maps 
with multiple clumps of matter and we tested the effects of the 
smoothing on the ellipticity of the mass distribution. In prac¬ 
tice each clump was generated with a NFW prohle with c - 4 
and = 5.0 10^^ Mq (associated to = 60" = 360 kpc, 
^200 = 1.4 Mpc and M 200 = 3.9 lO''^ Mq). Three clumps were 
aligned along a line and regularly spaced, with separation rang¬ 
ing from 40", 60" and 80" (250 kpc to 500 kpc at z ~ 0.5) 
between the clumps. The 2D mass maps were then smoothed 
with 2 different kernels, with cr = 30" and cr — 60", corre¬ 
sponding respectively to a ICF of 75" and 150" in the mass re¬ 
construction. The ellipticity of the simulated mass distributions 
was measured with the same method as for the clusters maps, 
with the centering fixed on the central mass peak, for both the 
smoothed and unsmoothed distributions (Fig|^. The relevant el- 
lipticiy was measured in radii between 100" and 150" for the 
observed clusters. The measures on the simulated clusters show 
that the apparent ellipticity is typically decreased by a factor 2 
when applying this severe smoothing. We did not attempt to cor¬ 
rect more accurately the measured ellipticities because the ellip¬ 
ticity attenuation should also depend on the mass profile, which 
is not constrained enough in this study. A higher background 
galaxy density would have allowed sharper mass reconstructions 
but this is out of reach with ground-based wide held imaging and 
requires data on the quality of HST 
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Table 4. General properties of the weak-lensing mass maps for the cluster sample: signal-to-noise ratio of the central peak of the 
mass map, ellipticity and position angle (PA), projected distances of the mass density peak from the brightest cluster galaxy (BCG) 
and the peak of the X-ray emission, and distance between the peak of the X-ray emission and the BCG. The elliptical parameters are 
ht from a mass isocontour drawn at the 3cr level, except for 3 clusters (see text for details). Position angles (PA) are given in degrees 
counter-clockwise with respect to the NS axis. Shifts between the BCG and the mass peak smaller than 30" can be considered as 
non signihcant and values larger than this limit are marked in bold. Note that in most cases, the position of the BCG does not exceed 
6" (or about 40 kpc) from the X-ray peak. 


Cluster 

Peak S/N 

Ellipticity 

6 

PA 

(deg.) 

Distance 
(Peak, BCG) 

Distance 
(Peak, X-rays) 

Distance 
(X-rays, BCG) 

MS 0015.9+1609 

6.5 

0.26 + 0.11 

+95 + 9 

43" 

49" 

6" 

(38 kpc) 

MS 0451.6-0305 

4.3 

0.38 ± 0.09 

+ 134 + 6 

33" 

27" 

3" 

(19 kpc) 

RXCJ0856.1+3756 

5.5 

0.20 ± 0.07 

+54 + 9 

6" 

4" 

4" 

(22 kpc) 

RX J0943.0+4659 (*) 

6.4 

0.18 + 0.05 

+88+12 

72" 

7" 

75" 

(407 kpc) 

RXC J1003.0+3254 (*) 

5.5 

0.35 + 0.03 

+21+9 

82" 

94" 

151" 

(830 kpc) 

RXJl 120.1+4318 

2.7 

0.15 + 0.05 

+42 + 13 

69" 

63" 

28" 

(190 kpc) 

RXC J1206.2-0848 

7.0 

0.25 + 0.04 

+74 + 4 

5" 

3" 

2" 

(11 kpc) 

MS 1241.5+1710 

4.4 

0.16 + 0.06 

+ 149 + 21 

13" 

16" 

1" 

(6 kpc) 

RX J1347.5-1145 

9.8 

0.13 + 0.05 

-19 + 20 

4" 

< 1" 

< 1" 

(< 6 kpc) 

MS 1621.5+2640 

5.4 

0.27 + 0.06 

+ 149+13 

26" 

42" 

8" 

(45 kpc) 

RX J2228.5+2036 (*) 

4.7 

0.40 + 0.07 

+ 104 + 5 

40" 

56" 

4" 

(22 kpc) 


(*) double cluster or merger 



Fig. 5. Ellipticity of simulated clusters formed by 3 mass clumps 
with NEW profile, and smoothed with 2 different kernels. The 
continuum line is the ellipticity measured on the unsmoothed 
data. Eor comparison the ellipticity of the observed clusters was 
measured in radii ranging from 100 to 150" approximately. 


The results of the elliptical fitting of the cluster mass dis¬ 
tribution are presented in Table Half of the clusters has a 
low ellipticity (5 clusters with edm < 0.2) while the other 
half is clearly elliptical (6 clusters with sqm > 0.2). The el¬ 
lipticity distribution of the sample shows an weighted mean of 
(e) = 0.25 + 0.12. It appears narrower than that of other sam¬ 
ples like LOCUSS ( |Oguri et al.||2010 i. But we have a smaller 
number of clusters and the difficulties to provide weak lensing 
maps at redshift 0.5 are stronger. No value exceeds 0.40, con¬ 
trary to what is expected for such a cluster sample and none of 
the clusters is really circular (edm < 0.1) in their extended re¬ 
gions. The smoothing process, needed to get an acceptable mass 
map at high redshift, is certainly the cause of the lack of high el¬ 
lipticity clusters in our sample (Eig|^. In conclusion, even if we 
can not draw firm conclusions on the ellipticity distribution of 


the dark matter from our sample, we show that we remain com¬ 
patible with standard expectations and previous works. A similar 
attempt to test the evolution of the ellipticity of the X-ray gas dis¬ 
tribution was proposed by Maughan et al. ( |2008[ ) with Chandra 
data. They did not hnd any change in the ellipticity distribution 
of the high-z sample (> 0.5) compared to their low-z sample, 
contrary to other morphological parameters like the slope of the 
surface brightness prohle at large radius. 

We also compared the distance between the main mass peak 
and the location of the BCG to check the consistency between 
the positions of the dark matter peak and the light peak. However 
the uncertainty in the centering of the mass distribution is high 
and strongly depends on the S/N ratio of the main mass peak. 
As demonstrated by Dietrich et al. ( |201 l| l with numerical simu¬ 
lations, the shape noise in the weak lensing map reconstraction 
combined with the smoothing process generate an offset distri¬ 
bution with a mode as large as 0.3' and median values up to 1' for 
typical ground-based observations. In our case, the clusters with 
the best map reconstruction and the highest S/N in the central 
peak (> 5) are also those for which the position of the mass peak 
matches the position of the BCG as well as the centroid of the 
gas distribution (Eig. |A.l| to A. 11 in the Appendix). This is gen¬ 
erally valid, except for the clusters RX J0943.0H-4659 (merging 
cluster) and MS 0015.9H-1609. This last cluster deserves some 
comment because previous weak lensing modeling of the cen¬ 
tral area, using HST/ACS images, point towards a mass centre 
well centered on the three brightest galaxies ( |Zitrin et al.|201 l| l. 
Our mass reconshTiction suffers from the proximity with a bright 
star and its halo which distorts the shape distribution of the faint 
galaxies. 

More interesting is the distribution of the X-ray/BCG off¬ 
set. It clearly appears bi-modal with most of the clusters hav¬ 
ing an offset smaller than 50 kpc, and three outliers. Among 
these three clusters, two are merger systems and the last one, 
RX J1120.1 h- 43 18 is rather poorly defined in its centre. This is a 
similar trend as found by Sanderson et al.|p009 1 in the LOCUSS 
sample of low redshift clusters. This X-ray/BCG offset is a good 
indicator of the dynamical state of the cluster and is highly cor- 
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Fig. 6. Ellipticity of the optical light distribution versus the el- 
lipticity of the dark matter. The vertical line separates the cluster 
sample between the “circular” clusters with e < 0.2 and the “el¬ 
liptical” or the irregular ones. 


related with the strength of the cooling core in the centre. The 
comparison between the X-ray and the mass peaks is more un¬ 
certain due to the limitations mentioned above. But it shows a 
similar trend, at least for the clusters with a mass peak detected 
with a high enough significance. 


5. Mass and light distributions 


5.1. Comparison between light and mass 2D distributions 

Fig. 1^ shows the correlation between the ellipticity of the dark 
matter and that of the stellar light distributions. As for the com¬ 
parison between the cluster light and the BCG ellipticities, we 
measured a weighted mean < sqm - sught >= -0.1 ± 0.02 and 
the mean orthogonal deviation to the sdm - flight line is 0.08 for 
the 11 clusters. Dispersion is slightly higher than for the light- 
BCG comparison but again we find a good concordance between 
both ellipticities, and a tendency towards a better agreement for 
elliptical clusters than for circular ones. We remind that these 
measurements correspond to large scale morphologies so they 
are more sensitive to substructures which can be found at the 
Mpc scale. Weak lensing morphology can also be disturbed by 
additional mass halos projected on the line of sight but not phys¬ 
ically related to the clusters. Similar conclusions were reached 
by Oguri et ar| ( 2010| l from a very similar study. 

To study the possible alignment effect between the light and 
the dark matter distributions, we represented the position angle 
of the optical light versus the position angle of the dark matter 
distribution (Fig. |^. To better visualize the shift with respect to 
the y - X line, we also computed the distance between the data 
points associated to the clusters and the 1;1 line in that plane. 
The position angles of both distributions are strongly correlated; 
the mean difference < APA > - < PA(DM) - PA(light) > = 
h- 3 ± 3° degrees and the average deviation to the 1:1 line is 28°. 
In addition, as shown in Fig. these differences tend to vanish 
when the ellipticity increases. We suspect that large differences 
at small ellipticities are partly due to biases in the processing of 
the elliptical fits or to the influence of further substructures at 
large radius. However |Oguri et¥r] ( |2010] l show in their detailed 
study that elliptical fits of weak lensing maps are robust when 
they use similar radii of 400 to 800 kpc. So we are confident 
that in the present study the “light traces mass” assumption is 


Fig. 7. Orientation of the optical light distribution versus the 
dark matter. The “circular” clusters with an ellipticity e < 0.2 
are marked with a cross. 
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Fig. 8. Orthogonal deviation (in degrees) from the 1 : 1 line for 
the dark matter orientation (PA(DM)) compared to the light dis¬ 
tribution orientation (PA(light)) versus the ellipticity of the dark 
matter distribution. The more elliptical this distribution is, the 
better the alignment between light and dark matter is. 


valid when clusters are in quiescent phases of their evolution. 
Departures from this assumption occur when the clusters enter 
merging processes and when interactions between large clumps 
of matter globally perturb their dynamical equilibrium. 

5.2. Cluster mass-to-light ratio 

The mass-to-light ratio is a quantity that has been widely studied 
at every mass scale, from single galaxies up to rich and massive 
galaxy clusters. Its variation across the mass range allows, for 
instance, to highlight physical processes that affect the baryonic 
component of massive structures, like star formation or galaxy- 
galaxy interactions in large dark matter halos ([Carlberg et al. 
|1996[|Marinoni & Hudson|2002t [Giodini et al.|2009| l. In Paper 

1, we analyzed the correlation between mass and luminosity and 
found a logarithmic slope of 0.95 + 0.37 for the mass-luminosity 
scaling law, compatible with a constant M/F ratio. In the present 
paper, we used 2D projected mass instead of the 3D mass 
used in Paper 1. Both masses differ mostly by a scale factor 
because all clusters are assumed to have the same concentra- 
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Fig. 9. Mass-to-light ratio versus the total luminosity for the 11 clusters of the sample (color points) and the 11 clusters added 
from the Bardeau et al. (2007) sample. The average value MjL - 166 in solar units is drawn as the solid line. 


tion parameter. This projected mass is the right quantity to be 
compared to the total projected luminosity. So we correlated 
the 2D mass with the total luminosity of early-type galaxies, 
computed inside the radius /?200 and corrected from the non- 
detected part of the luminosity function. We obtain an average 
ratio {MjL) — 160 + 60 h-jo (MIL)q, with values ranging from ~ 
68 to 235 (Table|^and Fig.|^. Our results are in excellent agree¬ 
ment with the values obtained by Bardeau et al. (2007) who used 
a similar methodology to derive the weak-lensing masses and the 
optical luminosities for a sample of clusters at lower redshifts 
{z ~ 0.2). They found an average ratio of {MjL) - 170 + 67 
for the same quantities, and their results are very similar to the 
present ones. Because the methodology they used is so close to 
the present one, we merged the two samples although the clus¬ 
ters differ in redshift. From the 22 clusters we hnd an average 
ratio {MjL) — 166 + 62 h-jQ (MfL)Q. Comparisons with other 
samples are made rather difficult due to the several methodolo¬ 
gies used to derive optical luminosities ( Popesso et al.|2007] l and 
cannot be discussed further. 

In addition, although the mass interval of the sample is rather 
limited, ranging from 6 x lO''^ Mq to 2.5 10'^ Mq, we tried to 
consider the possible variation of the M/L ratio with mass. This 
is rather speculative and limited by the fact that the clusters with 
the lowest mass are also those with the largest uncertainty in 
the weak lensing peak detection and that the 2 quantities are 


correlated. Following the same procedure as described in Foex 
et al. (2012|l, we htted the M versus L relation using a linear re 


gression in the log-log plan with the orthogonal BCES method 
which takes into account errors in both directions and provides 
a statistical dispersion around the ht crs,at as well as the intrin¬ 
sic dispersion cr - int. We did not hnd any signihcant depar¬ 
ture from a constant ratio between the 2 quantities, with a slope 
a - 0.945 ± 0.31, (Tint — 0.14 and cTstat — 0.11 for the disper¬ 
sions of the ht in the log-log space. With the present data we 


do not hnd any departure from a constant M/L ratio independent 
of the total cluster mass, although previous works have found 
a signihcant increase of the M/L ratio for massive clusters of 
galaxies compared to lower mass clusters and groups (|Popesso 


|et al.|20()7{|Andreon|2010[ ). But this was obtained with samples 
spanning a much larger mass interval than the present one. The 
physical origin of this situation is still controversial but it is usu¬ 
ally understood in terms of a decrease of the star formation ef- 
hciency with increasing halo mass (|Springel & Hernquist|2003 
|Lin et al.|2003l ). 


6. Summary and conclusions 

The cluster sample presented in this paper is limited to 11 clus¬ 
ters only, at a redshift z ~ 0.5. They were selected by their X-ray 
emission and are part of the representative sample EXCPRES 
which includes 20 clusters in the redshift range 0.4 < z < 0.6. 
Because of additional criteria used to optimize the weak lensing 
detection and analysis, the present sample is not any more rep¬ 
resentative of the cluster population at intermediate redshift, but 
it forms a sub-sample with the brightest X-ray luminosity. We 
summarize below the properties of the clusters which have been 
explored in the paper: 

- We provide for each cluster the total luminosity of the clus¬ 
ter, after a careful identihcation of cluster members, and sev¬ 
eral morphological parameters of the light distribution. We 
hnd good correlations between the ellipticity of the BCG and 
the global light distribution in terms of ellipticity and orien¬ 
tation. But whether the BCG has a bright and extended enve¬ 
lope of cD-type or not does not make signihcant differences 
in the general optical properties of the clusters. 

- The weak lensing mass reconstruction is done for each clus¬ 
ter, although the peak detection is at low signihcance in a 
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few cases (one cluster detected at less than 3cr). The aver¬ 
age ellipticity of the mass maps is < e >= 0.25, a value 
compatible with similar estimates at lower redshift. No evo¬ 
lution of the average ellipticity of the clusters or the fraction 
of high ellipticity mass distributions is detectable in our data. 
We also explore the distance between the mass peak, the lo¬ 
cation of the BCG and the X-ray centre for each cluster. The 
position of the mass peak is the most uncertain and is lim¬ 
ited intrinsically by the low density of background galaxies 
in the mass reconstraction. On the contrary, we hnd good 
agreement between the location of the BCG and the X-ray 
centre, especially for regular clusters. As expected the most 
discrepant clusters are those with the most disturbed mor¬ 
phology or clear signs of dynamical perturbations. 

- The mass-to-light ratio distribution shows excellent agree¬ 
ment with previous measures done with a similar approach 
and the average M/L ratio is found; {MIL) - 160 + 
60 h-jo {MIL)q. Previous studies were done at lower redshift 
and we do not hnd signihcant sign of evolution, as expected 
for this intermediate redshift bin. 


These properties point towards a general picture of the clus¬ 
ters for which “the light follows mass” paradigm is the main 
driver. This good coherence is valid both in the central parts of 
the clusters and at large scale, as demonstrated with the weak 
lensing mass reconstructions. Going into more details, we tried 
some attempts to separate the sample in two classes as was done 
previously for other cluster samples like LOCUSS or the CCCP 
( Smith et al.|2005[ Mahdavi et al.|2013] l. A majority of clusters 
are regular and follow the main correlations and 3 or 4 outliers 
are identihed as non-relaxed clusters or merger systems (namely 
RX J0943.0-1-4659, RXC J1003 .Oh- 3254, RX J2228.5-h2036 and 
possibly MS 1241.5 h- 1710). To better quantify departures from 
the assumption of regularity it will be important to study and 
better understand the inhuence of substructures and the role of 
triaxiality. It is a natural consequence of structure growth driven 
by self gravity of Gaussian density huctuations ([Limousin et al. 


2013|l but up to now it has been mostly neglected, for simplicity. 


We now have in hands a good understanding of the tracers of the 
distribution of the different components in clusters so it is timely 
and appropriate to address this issue in details. This would al¬ 
low to improve mass measurements and the understanding of 
the mass growth of structures as massive as clusters of galaxies. 
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Appendix A: Individual properties of clusters 

The 11 clusters presented in the sample are all bright X-ray clus¬ 
ters. Some of them also present specific optical properties or 
are already known as strong gravitational lenses. We review in 
this section the properties of the clusters, mostly in the optical. 
Their X-ray properties will be presented in a companion paper 
(Arnaud et al., in preparation). 


A.1. MS 0015.9+1609 iz^QMl) 

This very rich cluster has been studied for many years, since 
the identification of a high fraction of red galaxies in its pop¬ 
ulation ( Koo|| 198T| |. Included in the CNOC cluster sample, its 
spectroscopic survey was presented in [Ellingson et al.| (|1998|l 


with more than 180 objects observed spectroscopically. The re¬ 
sulting velocity dispersion crios - 1127) ;}^^ km s~^ is a high 
value consistent with its galaxy richness ( Borgani et al.||T999 i. 
MS 0015.9+1609 is one of the brightest and most distant X-ray 
cluster included in the EMSS sample ( |Gioia & Luppino|1994) l. 
It is also part of the highly luminous X-ray clusters identified 
in the MACS sample at redshift larger than 0.5 ([Ebeling et al.| 
|2007l ) and is identified as MACS J0018.5+1626. 


The weak lensing properties were described by Smail et al. 


(1995 1 and then by Clowe et al. (20001. The authors found 
a rather low signal and therefore a total mass not consistent 
with the optical velocity dispersion of the galaxies. More re¬ 
cently Hoekstra ( 2007| l re-analyzed a large sample of clusters 
observed in good seeing conditions at CFHT and found for MS 
0015.9+1609 a total mass described by a SIS with cr = 11642*^* 
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‘ M,. 


km s^' or by a NEW profile with M 200 = 27.0!® ° ivio 

Note that despite the high mass value of the cluster, no strong 
lensing features were detected in HST images (Sand et al. 2005 1 . 
More recently and thanks to a detailed analysis of HST/ACS im¬ 
ages, Zitrin et al. p011| identified 3 systems of multiple images, 
though not yet confirmed spectroscopically. They were used to 
provide a lensing model of the mass distribution in the centre of 
the cluster. 

There is no dominant central galaxy in this cluster but a chain 
of bright ellipticals, giving a significant elongation in the galaxy 
distribution. This elongation was confirmed in the weak lensing 
map provided by [Zitrin et al. (2011 1 on the central area of the 

lipticity of the mass distri- 
We we suspect that 


cluster. In our wide field map, the e 
bution does not clearly appear (Fig. [A.l 
the bright star which is close to the cluster centre prevents a cor¬ 
rect study of the cluster mass map obtained from weak lensing 
reconstruction. 

MS 0015.9+1609 is embedded in a large scale structure 
of the size of a supercluster, identified spectroscopically by 
Connolly et al. ( 1996|l. At least 3 clusters lie within less than 


30 Mpe form each other, and a long and massive filamentary 
structure crosses the cluster in the same direction as the galaxies 
elongation ( Tanaka et al.|2007 2009| l. The weak lensing recon¬ 
struction presented in the present paper is focused on the central 
area around the cluster only, but we checked that most of the 
structures spectroscopically identified by [Tanaka et al.[ ( 2007| l 
are also visible in our global mass map. This may be the case 
for the South-West elongation seen in the mass map displayed 
in Fig. A.l a Further work is in progress to better quantify these 
coiTelations. 


A.2. MS 0451.6-0305 (z = 0.537) 


This cluster is the most X-ray luminous cluster in the EMSS 
catalog ( [Gioia & Luppino[[1994[ ) and is also part of the CNOC 
sample. Intensive spectroscopic follow-up of the galaxies pro¬ 
vided more than 100 spectra of cluster members (Ellingson et al. 


TTT 


1998[) and a line-of-sight velocity dispersion of cr/oj = 1002_gj 


km s ' (Borgani et al. 19991. Weak lensing masses measured by 
Clowe et al. ( 2000[ ) are roughly compatible with this value as 
well as those obtained by [Hoekstra et al.[P012| l. Our own mea¬ 
surements are 50% higher, but they remain compatible within 
the uncertainties (|Foex et al.|20121l. The cluster is also identified 
as MACS J0454.1-0300. 

A few thin and elongated features were suspected as strong 
lensing candidates by Luppino et al. ( T999]l_^and spectroscop¬ 
ically confirmed later by Borys et al. (20041. Interestingly a 
SCUBA detection of an extended source in the cluster centre 
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lead to the identification of an ERO pair, triple imaged (|Chapman| 
|et al.|2002||Takata et al.|2003[|Berciano Alba et al.|2010[ ). These 

features point towards the bright central galaxy as the centre of 
the mass distribution. The observed elongation of the weak tens¬ 
ing mass reconstruction (Fig. |A.2| l is well correlated with the 
global elongation of the light distribution in the SE/NW direc¬ 
tion. This is also true for the orientation of the BCG. The latest 
strong lens model presented by Zitrin et ahlpOl 1[ ) indicates that 
the central mass distribution is highly elliptical with an orienta¬ 
tion that matches the SE/NW elongation of the cluster at large 
scale. 


|2007[ ). Nothing was really known on the optical properties of this 
cluster, which displays a bright galaxy in its centre and a rather 
loose distribution of cluster members. Another bright galaxy is 
located 2.3' South-West, with similar properties. It is centered on 
a secondary peak in the X-ray gas distribution and the mass map 
is centered in between these two galaxies. But the bi-modality 
of the cluster is more visible on the galaxy distribution than in 
the mass map (Fig. |A.5| l which is limited by its spatial resolu¬ 
tion. We suspect that this cluster results from the merging of two 
sub-clusters and all conclusions regarding RXC J1003 .Oh- 3254 
in the global analysis of the sample must be taken with caution. 


A.3. RXC J0856.1+3756 (z = 0.411) 

This cluster is part of the NORAS sample (Northern ROSAT 
all-sky galaxy cluster survey), a pure X-ray selected sample 
(Bohringer et al. 2000 1 . It was included in the EXCPRES 
sample because of its high X-ray luminosity. It is an opti¬ 
cal bright cluster i dentihed in the SDSS-DR6 ([WHL2009] 
J085612.7H-375615, ( |Wen et al.|2009| l), with a redshift measure¬ 
ment of the BCG at z = 0.411. The cluster displays a well de- 
hned and regular luminous over-density dominated by a bright 
and extended cD galaxy (Fig. |A.3| . The mass map also presents 
a very regular aspect around its centre and provides a coherent 
picture of a relaxed cluster. 


A.4. RX J0943.0+4659 (z = 0.407) 

This cluster is also known as Abell 851 or Cl 0939-1-4713. It 
is the only Abell cluster of our sample. High resolution HST 
images of the centre revealed a large population of blue galax¬ 


A.6. RX J1120.1+4318 {z^ 0.612) 


ies and a lot of merging galaxies (Dressier et al. 19941. Seitz 
et al. ( 1996| l used this deep HST/WFPC2 image to identify a 
few lensed objects but no highly magnihed gravitational arcs 
were detected. X-ray observations of A851, first with ROSAT 
and more recently with XMM- Newton, showed a very perturbed 
distribution with pronounced substructures and evidences for a 
dynamically young cluster ( |De Filippis et al. |2003| l. Tentative 
2D spectro-imaging lead to the identification of a hot region be¬ 
tween the 2 main sub-clusters, a characteristic of a major merger 
in an early phase. 

The galaxy distribution is complex, with a high galaxy den¬ 
sity in the central area. It can be separated into two clumps which 
trace the cluster interaction and are coherent with the gas dis¬ 
tribution. Several bright galaxies dominate the light distribution 
and are more concentrated in the South-West extension of the 
cluster (Fig. |A.4[ ). On the contrary the weak lensing mass map is 
surprisingly regular with only one main structure but elongated 
along the direction of the interaction. The separation between the 
two X-ray peaks is 50", well below the resolution of the mass 
map. So with the present data there is no chance to have a more 
detailed view of the mass distribution at a scale where the phys¬ 
ical processes of the cluster merger could be identified. Deeper 
imaging should be necessary to go further in this analysis. Due 
to the high evidence for merging processes, this cluster was re¬ 
moved later from the EXCPRES sample, but as optical data were 
obtained in good conditions, we kept it in the present sample. 


A.5. RXC J1003.0+3254 (z = 0.416) 

The cluster was identified initially by its X-ray extended emis¬ 
sion in the NORAS sample ( [Bohringer et al.| 20001 and it was 
later re-detected in the 400d ROSAT sample^Burenin et al.j 


This cluster belongs to the Bright SHARC survey ( jRomer et al. 
|2000| ) and it was included in the WARPS II catalog (the Wide 
Angle ROSAT Pointed Survey, [Horner et aL]p0081 l). The clus¬ 
ter was observed with XMM- Newton and analyzed by Arnaud 
jet al.j ( |2002| ) who found a regular X-ray emission with a spher¬ 
ical morphology. They also claim that no cooling flow or cen¬ 
tral gas concentration is present in this cluster, which is con¬ 
sistent with the cooling time being larger than the age of the 
universe at this redshift. Indeed with its redshift z = 0.612, RX 
J1120.1H-4318 is the most distant cluster of the EXCPRES sam¬ 
ple. The light distribution of cluster members shows a East-West 
elongation which was also been measured in the Chandra X-ray 
map ( Maughan et al.|2008| l. But the ellipticity is rather small and 
does not attenuate the regular morphology of the cluster which 
is clearly in a relaxed phase. Unfortunately the lensing signal in 
RX J1120.1-1-4318 is barely detected, at less than 3cr (Fig. |A.6| . 
So it is difficult to draw any conclusions on the mass distribution 
in the cluster. Even the shift between the mass peak and the light 
peak can not be considered as significant. 

A.7. RXC J1206.2-0848 (z = 0.441) 


This cluster is one of the brightest clusters of the REFLEX sam¬ 
ple (Bohringer et al. 2004| l. It belongs to the MACS sample 
( Ebeling et al.|2010 1 as MACS J 1206.2-0847 and is part of the 
CLASH sample ( [Postman et al.||201^ . It displays a bright and 
spectacular arc system, initially spectroscopically observed by 


Sand et al. ( [2004[ ) and confirmed more recently by Ebeling et al. 
( [2009[ l ata redshift z = 1.036. A detailed analysis of the central 
mass distribution was done both with strong lensing and X-ray 
data, giving a discrepancy of a factor 2 between the two mass es¬ 
timates. But the X-ray distribution of the gas shows some signs 
of merging processes in the centre which could explain this dis¬ 
crepancy. Similar trends have already been noticed in other clus¬ 
ters like A1689 for example ( P^imousin et al.|2007[ l. 

The weak lensing mass distribution is well peaked, with a 
regular shape and a central concentration which fits the lumi¬ 
nous mass as well as the X-ray mass (Fig. [A.7[ ). Note that the 
central galaxy is also a bright radio source with a steep spectrum 
( Ebeling et al.|2010j l. 

[Umetsu et al. ( 2012| l recently did a comprehensive analysis 
of this cluster combining weak and strong lensing derived from 
wide-held Subaru imaging and HST observations. Their mor¬ 
phological analysis of both the reconstructed mass map and light 
distribution reveals the presence of a large-scale structure around 
RXC J1206.2-0848. The orientation of this structure matches 
the position angle of the BCG and that of the cluster light distri¬ 
bution and projected mass map. The ellipticity they derive for 
the latter is somehow larger than ours, but we obtain consis- 
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tent results for the light distribution. The overall shape of RXC 
J1206.2-0848 indicates that light follow mass up to the large 
scales of the cosmic web. 


A.8. MS 1241.5+1710 (z = 0.549) 


part of the CNOC sample and was observed spectroscopically 
with m ore than 100 cluster redshifts available (|Ellingson et al. 


This is the brightest cluster of the REFLEX sample (Bohringer 
et al.| 2004l and is part of the CLASH sample ( [Postman et al. 
2012^. It presents a spectacular strong lensing system detected 
by [Schindler et ah ( 1995| l. It was also identified as a cluster 
with a strong central cooling flow ( Allen|2000l l, feeding a pow¬ 
erful radio source ( [Pointecouteau et al.|2001| l A detailed com¬ 
bined analysis of the strong and weak lensing effects ( [Bradac 
et al. 2005a|b| l lead to a very accurate view of the dynamical 
status of the cluster in the inner regions: the cluster presents a 
mass concentration centered on the BCG with some extension 
to the SW and many evidences of sub-cluster merging. But RX 
J1347.5-1145 is definitely not a major merger. After some con¬ 
troversy, the different mass estimates seemed to converge, espe¬ 
cially those measured close to the centre using the strong lensing 
features ( [Halkola et al.[[2008} [Bradac et al.[[2008| ). But a factor 
of 2 remains between the X-ray and the weak lensing masses 
at large radius ( [Fischer & Tyson[|l997[ Kling et ar|[20B5[ Gitti 


et al.|2007 1 . In this context, our mass map confirms the previous 
results and does not bring new evidences on the mass distribu¬ 
tion (Fig. [A.9[ l. It is mostly used in order to check and validate 
our weak lensing procedure before applying it to other clusters. 


Similar results were published on RX J1347.5-1145 by Hoekstra 
[et al.[ ( 20T^ with the same CFHT data. Fortunately they obtain 
very similar mass measures. 

Several studies ( [Lu et al.[2010[[Verdugo et ah 2012) revealed 
that RX J1347.5-1145 is embedded in a large-scale structure, ex¬ 
tending up to 20 Mpc in the NE-SW direction. Our luminosity 
map confirms the existence of several over-densities on a large 
scale, aligned along this direction. The main orientation of the 
cluster light distribution also follows the same direction. As for 
RXC J 1206.2-0848, this cluster supports the picture of the cos¬ 
mic web where massive clusters are fed by filaments whose ori¬ 
entation matches the global morphology of the central node. 


w 


1997(1. The velocity dispersion is not very high {crios - 839_53 


As part of the EMMS sample ( Gioia & Luppino[1994[ l, this clus¬ 
ter was also observed in optical, but no significant strong lensing 
feature was detected ( [Luppino et al.|1999] l. The luminosity dis¬ 
tribution is complex, with a southern extension possibly related 
to the main cluster. However the mass distribution does not show 
a similar trend nor the X-ray gas. Both are regular and centered 
in the BCG embedded in a bright and extended envelope. So firm 
conclusions are difficult to draw because of the low S/N ratio of 
the mass map (Fig. [A.8[ l. The second over-density of galaxies 
could also be due to some contamination along the line of sight. 
Deeper and multi-color images would be necessary to confirm 
the reality of an in-falling substructure on the main cluster. 

A.9. RX J1347.5-1145 (z = 0.451) 


imi/s, [Borgani et ali] ( [1999| l). More recently [Hoekstra[ ( (2007[ ) did 
a very accurate weak lensing analysis and his results are in good 
agreement with the dynamical mass estimate. It is also consistent 
with the X-ray mass obtained with ROSAT ( [Hicks et al.|2006[ l. 
With our weak lensing mass reconstraction, we And a mass dis¬ 
tribution rather elongated and coherent with the light distribu¬ 
tion. The large shift between the mass and light peaks is more 
probably an artifact than a real one (Fig. [A.10l l. 

A.11. RX J2228.5+2036 {z = 0.412) 


This cluster is part of the NORAS sample ( [Bohringer et al.[20()0 1 
and also belongs to the MACS sample (MACS J2228.5-1-2036, 
[Ebeling et al.[ ( [200'7] l). Because it is at low galactic latitude, very 
few optical observations are available. Our weak lensing recon¬ 
struction is rather uncertain (Fig. A.ll[ l and possibly flawed be¬ 
cause of the large number of bright stars in the held of view. 
However, in addition to its strong X-ray emission, this cluster 
was detected for its SZ signal, allowing one of the first combined 
analysis between the X-ray and the SZ signals ( [Pointecouteau 
[et al.[|2002j [Jia et al.][2008| l. Both confirm that the cluster is 
quite massive and dynamically perturbed. The weak lensing map 
shows a poor signal close to the cluster centre but suggests that 
the cluster has an elongated shape. This is also valid for the com¬ 
plex light distribution. The most convincing feature is a galaxy 
clump located in the South-West direction, detected on the mass 
map with higher significance than the main cluster. It is associ¬ 
ated with a galaxy excess centered on a bright elliptical galaxy 
with similar magnitude as the cluster BCG. We suspect that this 
clump is at a similar redshift as RX J2228.5H-2036 and may 
be the cause of a future major merger with RX J2228.5-1-2036. 
Surprisingly there is no X-ray counter-part to this clump. 


A.10. MS 1621.5+2640 (z = 0.426) 

As part of the EMSS cluster sample (Gioia & Luppino 1994| l 
the cluster was rapidly identified as a strong lens with a nice 
gravitational arc located around a radio galaxy which is not the 
brightest cluster galaxy ( Luppino et al.|1999] l. No spectroscopic 
redshift is presently available for the arc, although its lensed na¬ 
ture is not in doubt ( [Sand et al.|2005[ ). MS 1621.5H-2640 is also 
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Right ascension 



Fig.A.l. Up: 15' X 15' insert of the cluster field extracted from the full r' MegaCam image. The thick red contours show the 
mass distribution derived from the 2D weak lensing analysis. The contour levels are linearly spaced in cr of the mass reconstruction, 
starting at 2cr. The thin blue contours come from the X-ray map obtained with XMM- Newton. The image was hltered with wavelets 
and the contours are scaled logarithmically. Bottom left: same mass isocontours overlaid on the galaxy luminosity distribution where 
cluster members are selected within the cluster red sequence and iriy < 23. Bottom right: true color image of the cluster centre from 
g', r', i' combination. The field of view is 3' x 3'. 
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Fig. A.2. Same as Fig. 


A.l 
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Fig. A.4. Same as Fig. 


A.l 
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Fig.A.5. Same as Fig. 
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Fig. A.7. Same as Fig. 


A.l 
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Fig. A.9. Same as Fig. 


A.l 
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Fig.A.lO. Same as Fig. 
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